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1. Introduction 

A large part of the earth is covered by thin semi-transparent cirrus 
cloud (Wylie et al . , 1989). The cirrus results from the natural injection of 
moisture into the upper troposphere by deep convection (i.e., anvils) and from 
man-made moisture injected into the upper troposphere by jet aircraft. 

Although most cirrus cloud is se mi- transparent to infrared wavelengths, their 
heights, thicknesses, and spectral absorption properties must be known in 
order to retrieve atmospheric temperature and moisture profiles from the data. 

In this paper, an algorithm is developed for accounting for the radiative 
properties of semi- transparent cloud in the retrieval of vertical temperature 
and moisture profiles. The algorithm is to be applied to the NASA ER2 HIS 
data collected during the FIRE cirrus field program. The results of its 
application will be reported at the second FIRE Annual Meeting to be held July 
10-14, 1989 in Monterey, California. 

2. The Cloud Algorithm 

The retrieval of temperature and moisture profiles from spectral 
measurements of upwelling radiances are most efficiently achieved through the 
simultaneous solution of the radiative transfer equation of the form (Smith et 
al. , 1989) 
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where I is spectral radiance, v is wavenumber, T is temperature, f is the 
derivative of Planck radiance with respect to temperature computed for a 
"guess" temperature profile condition, r° is the total atmospheric 
transmittance (i.e., the transmittance due to all absorbing constituents) for 
the "guess" atmospheric conditions, r? is the transmittance due to 
individual absorbing constituents (e.g., CO 2 , ^0, and cloud). The subscript 
"s" denotes the surface value, N the number of individual absorbers to be 
considered, and the transmittance function is understood to pertain to the 
atmosphere between the instrument level (p*) and the level of interest (p) . 

The symbol 6 denotes a perturbation from the "guess" condition. For a 
constituent whose concentration is known apriori, <5T^ represents the deviation 
of the actual temperature profile from the "guess" profile. For all other 
constituents , 
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where ST - T-T°, is the total path length of the gas between the instrument 
and the level of interest. The significance of Eq. (1) is that it is linear, 
the only unknowns are the temperature perturbations, ST^p) and they can be 
determined by the linear inverse of (1) . 

Neglecting scattering, clouds can be treated in the very same manner as a 
molecular absorber, provided its spectral dependence is known. In this case, 
the total transmittance r°(i/j) is given by 
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where the "guess” cloud transmittance (e.g., i-3) has the form 


T cd^j) - r °3^j> - EXP ( - Q f ^j- r °) ) 

Pb'Pt 


(3) 


In (3), a is a constant dependent upon the total ice or water content of the 
cloud, f(i/,r°) is a spectral function dependent on the effective radii, r° , of 
the absorbing particles (ice or liquid water drops), P t and are the top and 
base pressures of the cloud. Assuming that the cloud is purely absorbing, a 
good approximation for f(u,r°) can be obtained by the Modified Anomalous 
Diffraction Theory (MADT) presented by Ackerman and Stephens (1987). In this 
case , 
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where Q ABS is the absorption efficiency, * is the size parameter 2nr°/X, where 
r is the effective particle radius and A is the wavelength, m is the index of 
refraction, and ^ is the complex part of the index of refraction. Figure 1 
shows a comparison between the absorption efficiency computed for ice clouds 
using modified anomalous diffraction theory and Mie theory. It can be seen 
that the difference in spectral dependence for large and small particles is 
explained with MADT. In fact, if the effective radius is a free parameter, as 
it is in the retrieval problem, an effective radius for Eq. (4) can be chosen 
which provides an even closer fit to the Mie calculations. 
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Thus, given (3) and (4), one has the following four cloud parameters, a, 
r°, P t , and . If we consider the cloud term of Eq. (1), we can write 
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where /3 = a/(P^-P t ) . Thus, given an initial guess of a, r°, P t , and P^ from 
which r°(i/^)can be specified, one could solve for the profile 
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As with the molecular absorber 
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Since 
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the cloud temperature profile, T C(i , 
Eq. (1), Eq. (6) can be solved for 6/3 
and Sp t . Given 6/3 and 6 p t , a new base pressure can be defined from 


T cd - T - - - [ 6P (— — ) - «P t ] • 

3p 0 ° 

If the true temperature profile, T, and 
are known from the increase solution of 


a a 

Pb “ Pt + — ( Pt° + + 

0 


I 8°+60 


(7) 


The solution for the absorber temperature profiles and cloud parameters depend 
upon an initial estimate of the particle radius r° and the particle 
concentration variable a 0 . It can be shown that for an isothermal cloud, and 
at a wavenumber void of molecular absorption (i.e., a "window") 
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In order to solve (8) for r° and a 0 , the initial cloud top pressure can 
be defined using the C0 2 slicing technique (Menzel, 1986; Smith and Frey, 
1988). Given the product f(i/ 0 ,r°)a 0 for two window wavelengths whose cirrus 
optical properties are expected to differ (e.g. , 8.5 and 11.5pm as shown in 
Fig. 1), r° can be defined from the ratio of Eq . (8) applied to the two 
wavelengths since a 0 cancels. Once r° is known, a° can then be obtained from 
Eq . (8) applied to either one of the two window wavelengths. Given a° an 

estimate of fi° can be obtained from an initial guess of cloud-base altitude 
p b ° ( i . e . , p° — a°/(p b °-p t °)). 

In summary, the following steps can be taken to solve simultaneously for 
temperature and water vapor profiles and cloud parameters from cirrus cloud 
contaminated radiance spectra: 

(1) Using the C0 2 slicing technique with the guess temperature profile, 
solve for an initial cloud top pressure, p° . 

(2) Using the 8 . 5pm and 11.5pm window radiances and the initial guess 
surface and cloud temperature conditions, use Eq . (8) to determine the ratio 
f(8.5pm, r°)/f (11 . 5pm, r°) . Use (4) to determine r° and then Eq . (8) applied 
to one of the two window wavelengths to determine a°. This step yields the 
initial guess cloud parameters a , /3° , p b ° . 

(3) Specify the initial guess cloud transmission profile using Eq . (3) 

where f (i/j , r°) is calculated using Eq , (4). 

(4) Solve for the profiles through the linear matrix inverse 

solution of Eq . (1). 

(5) Use Eq. (6) to solve for S/3 and 6p t . 

(6) Use Eq. (7) to provide p b . 

The entire process could be iterated until the cloud parameters cease to 
change from one iteration to the next. 

3 . Summary 

A method for direct and simultaneous solution for temperature and 
moisture profiles and cloud parameters is outlined. The method takes into 
account the spectral dependence of cloud emiss ivity/transmissivity which has 
been found to be a function of particle size. The results of the application 
of this technique to cirrus radiance spectra achieved with the HIS 
interferometer flown on the NASA ER2 during the FIRE will be reported at the 
meeting . 
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